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Abstract 

Ritt-Wu's algorithm of characteristic sets is the most representative for triangularizing sets of 
multivariate polynomials. Pseudo-division is the main operation used in this algorithm. In this 
paper we present a new algorithmic scheme for computing generalized characteristic sets by 
introducing other admissible reductions than pseudo-division. A concrete subalgorithm is de- 
signed to triangularize polynomial sets using selected admissible reductions and several effective 
elimination strategies and to replace the algorithm of basic sets (used in Ritt-Wu's algorithm). 
The proposed algorithm has been implemented and experimental results show that it performs 
better than Ritt-Wu's algorithm in terms of computing time and simplicity of output for a 
number of non-trivial test examples. 
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1. Introduction 



For solving systems of multivariate polynomial equations there are mainly three clim- 
ination approaches based o n r esultants, tria ngular sets, and Grobner bases (see, e.g., 
Kapur and Lakshman 19921 and IWarigll2001bl ). One of the b est known c oncepts for trian- 
gular sets is characteristic set, which was introduced first by iRitt (Il950h for prime ideals. 
Since 1980, W.-T. Wu has considerably developed Ritt's theory and method of charac- 
teristic sets by removing irreducibility requirements and designing efficient algorithms 
for zero decomposition of arbitrary polynomial systems. Ritt-Wu's method has been im- 
proved and extended by a number of researchers a nd has been succes sfully applied to 
many problems in science and engineering (see, e.g.. |Wu and Gaoll2007r) . 

To speed up the computation of characteristic sets JChouT i 19881 ) and later Chou and Gaol 
( 1990l ) introduced the notions of W-p rem and W-charac t eristic set. Characteristic sets 
and their complexity were studied by iGallo and Mishral (|l991allbh and new algorithms 
for computing characteristic sets with si mple exponen tial sequential and polynomial par- 
allel time complexities were presented. IWand (|l992) proposed an effectiv e strategy to 
improve Ritt-Wu's algorithm of characteristic sets. Gao and Choul ( 1992 ) presented a 
method based on Ritt-Wu's algorithm to deal with parametric polynomial systems. Cer- 
tain properties about a scending chains were s tudied and used to enhance the efficiency of 
Ritt-Wu's algorithm in I Gao and Chou ( 1993 ) . A complete i mplem e ntation of Ritt-Wu's 
method in the Maple system was reported in IWand (Il995h . IWand (l2001ah also general- 
ized Ritt-Wu's algorithm by means of one-step ps e udo-r eduction with strategies for the 
selection of reductends and optimal reductors. Q ( 20061 ) described a modified Ritt-Wu 
algorithm that can avoid redundant decompositions. 

Besides Ritt-Wu's, there are many ot h er effi cient methods for decomposing systems 
of multivariate polynomials. iKalkbrenerl (119931 ) introduced the notion of regular chain 
and presented a method for decomposing any algebraic variety into unmixcd-dimcnsional 
co mponents represented by regular chains. Regular chain was also defi ned independe ntly 
by lYang and Zhand (|1994l ) under the name of normal ascending set. iLazardl (fl99ll) in- 
troduced the concept of normalized triangular set (which is a special regular chain) 
and provided a method for uniquely decomposing the zero set of any polynomial set 
into regular zero sets of nor malized triangular sets. An o ther m ethod of triangular de- 



composition was proposed by Wang ( 19931 ) . Aubrv et al. ( 19991 ) studied the relationship 



between several different notions of triangular sets used in these methods and proved the 
equivalence of four notions (including regular chain). Comprehensive investigations o n 
triangular sets and trian gularizing alg o rithms were also carried out in iHubertl (l2003al) 
The r e ader m a y refer to Szantol Jl999 ); Aubrv and Moreno Maza (1999); Moreno Maza 
(|2000h : IWand |2000l l2001bl ): IChen and Moreno Mazal (|201ll ) and references therein for 
other related work on triangular decomposition of polynomial systems. Extensions of the 
theories and me thods in the algebra i c case to the differen t ial, difference, and other cases 
may b e found in fBoulier et all (ll995l):lLi and Wand (ll999h:lHubertl (l2000h:IChen and Gaol 
(|2003l ): lHubertl (|2003bl ) ; iBoulier et all (|2009l ); lGao et al.l (|2009l ) ; iBoulier et all (|2010l ) and 
references given therein. 

Pseudo-division is the main operation used in most of the above-mentioned methods. 
In the process of computing a characteristic set, one needs to construct polynomial re- 
mainder sequences (PRSs). The coefficients of intermediate polynomials in a PRS may 
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swell dramatically. As pointed out bv lCollins ( 19671 ) . the coemcients of polynomials in an 
Euclidean PRS (computed u sing pseudo-d ivision) may grow exponentially. An exponen- 
tial upper bound is given in iKnuthl (|l997l (27) in 4.6.1). There are several other PRSs, 
such as subrcsultant PRS in which the coefficients of polynomi als c an be determ ined 
by using certain matrices (see lvonzur Gathen and Lucking! 2003 and lMishralll993 ) and 
swell much more slowly than those in an Euclidean PRS. |lj dl98g) was the first who 
studied superfluous factors appearin g^ in the compu tation of characteristic sets. 

In this paper we follow the work of lWang (2001a) to present a new algorithmic scheme 
for computing generalized characteristic sets efficiently. To compute characteristic sets, 
one can replace pseudo-division by one-step pseudo-reduction, but this does not enhance 
the efficiency very much because extraneous fac tors for the pseud o-rests may be created 
(see the analysis in lSasaki and Furukawalll984l and IWang||2001al ). Our intention here is 
to design a reduction mechanism that can take advantage of the structure and properties 
of the polynomials under consideration and thus provide more flexibility for reduction 
strategics. A polynomial set which generates the same ideal as the input polynomial set is 
maintained, leading to the concept of generalized characteristic set. This set may contain 
polynomials of degrees smaller than the degrees of those in the input set and therefore 
may take less computing time for pseudo-reduction to 0. By introducing admissible re- 
ductions other than pseudo-division, it is possible to control the swell of coefficients of 
intermediate polynomials and to compute with smaller polynomials. Several strategies 
for concrete admissible reductions are adopted in a sample algorithm of the new scheme. 
Two versions of the sample algorithm have been implemented and compared with three 
algorithms implemented in Maple: two in the Epsilon package and one in the Groebner 
package. More details about the outputs of these algorithms for some examples are given 
in the appendix. 

The paper is structured as follows. Section 2 consists mainly of terminologies and no- 
tations about characteristic sets. In Section 3, we describe the main algorithmic scheme 
for computing generalized characteristic sets. A concrete subalgorithm is designed to tri- 
angularize polynomial sets by means of admissible reductions and effective elimination 
strategies and to replace the algorithm of basic sets (used in Ritt-Wu's algorithm). The 
termination and correctness of these algorithms are proved. In Section 4, discussions 
on concrete admissible reductions are given. Section 5 provides details for the sample 
subalgorithm. Section 6 contains an illustrative example and some experimental results, 
showing that our proposed algorithm performs better than Ritt-Wu's in terms of effi- 
ciency and simplicity of output for a number of non-trivial test examples. 



2. Preliminaries 

Let K. be any field. The notation K,\x] or K\x\, . . . ,x n ) stands for the ring of poly- 
nomials over K, in the variables x\, . . . ,x n . In what follows, we always assume that the 
variables are ordered as si < • • • < x n . 

Let F be an arbitrary non-zero polynomial in IC[x\. We use deg(F,Xk) and \c(F,Xk) 
to denote the degree and the leading coefficient of F with respect to (w.r.t.) the variable 
Xk respectively. The biggest variable that effectively appears in F is called the leading 
variable of F and denoted by lv(F). Moveover, ini(F) := lc(F,W(F)) and ldeg(F) := 
deg(-F, lv(-F)) are called the initial and the leading degree of F respectively. For any non- 
constant polynomial F, the index c of its leading variable x c is called the class of F and 
denoted by cls(F). The class of any non-zero constant is set to be 0. 



3 



Definition 1. An ordered polynomial set pi,. . . ,T r ] C K[x] \ K, is called a triangular 
set in JC[x] if lv(Ti) < • • • < lv(T r ). 

Definition 2. Let P and Q be two polynomials in K[x\. We say that P is reduced w.r.t. 
Q if deg(P, lv(Q)) < ldeg(<5); otherwise, P is reducible w.r.t. Q. 

Furthermore, let P = [Pi, . . . ,T r ] be a triangular set in K\x\. F is reduced w.r.t. P 
if for every i = 1, . . . , r, P is reduced w.r.t. Ti; otherwise, P is reducible w.r.t. P. The 
triangular set P is called an auto-reduced (or initial-reduced) set if for every i = 2, . . . , r, 
Ti (or ini(Tj)) is reduced w.r.t. [Pi, . . . , P;-i]. 

Definition 3. Let P and Q be non-zero polynomials in JC[x\. P is said to have a lower 
ranking than Q, denoted as P -< Q, if cls(P) < cls(Q), or cls(P) = cls(Q) and ldeg(P) < 
ldcg(Q). In this case, we also say that Q has a higher ranking than P and denote it as 
QyP. 

If neither P -< Q nor Q y P holds, then we say that P and Q have the same ranking 
and write P ~ Q. Moreover, P >- Q or P ~ Q is denoted as P £3 Q, and P ^ Q can be 
similarly defined. 

Definition 4. Let P = [Pi, ... , T r ] and S — [Si, ... , S s ) be two triangular sets in JC[x], 
P is said to have a lower ranking than S, denoted as P -< S, if 

(a) there exists an i < min(r, s) such that Pi ~ S\, . . . , Tj_i ~ Si-i and Ti ~< Si; or 

(b) r > s and Ti ~ Si, ■ ■ ■ , T s ~ S s - 

In this case, we also say that S has a higher ranking than P and denote it as S y P. 

If neither T y S nor P -< S, then P and S are said to have the same ranking and 
denoted as P ~ 5. In this case, one can easily know that r — s and Pi ~ Si, ■ ■ ■ , T r ~ S r . 

The order ^3 defined above is a partial order and can be used to order triangular 
sets. The following proposition indicates that any non-empty set of triangular sets has a 
minimal element. 



Proposition 5 (jWul 1 1986( 1. For any triangular set sequence Ti ^ pj ^ 7i ^ • • • , there 



exists an integer m such that for any i > m, Ti^Tm- 

An (initial-reduced) auto-reduced set in K.[x] is also called a (weak) non- contradictory 
ascending set, while [a] (a 6 /C\{0}) is called a (weak) contradictory ascending set. Both 
of them are jointly called (weak) ascending sets. 

Let V be a non-empty polynomial set in K[x\ and Q, denote the set of all the ascending 
sets that are contained in V; ri is not empty because the set of any single polynomial in 
V is an ascending set by definition. By Proposition 5, 51 has a minimal element, which is 
ca lled a basic set of P. This set can be computed, e.g., by the algorithm BasSet described 
in I Wand (l2001bh . 



Definition 6. Let P be a non-empty polynomial set K\x\. An ascending set M. is called 
a medial set of P, if M. C (P) and A4 has ranking not higher than the ranking of any 
basic set of P. 



Proposition 7 (|Wandl200l"bL Let P be a non-empty polynomial set in K,[x], M. = 



[Mi,...,M r ] a medial set of V with Mi $ K., and M a non-zero polynomial reduced 
w.r.t. M. Then for any medial set M' of V' = V U M U {M}, we have M' -< M. 
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For any f,G £ K\x] with G ^ 0, let prem( F, G, Xk ) and pquo(F, G, Xk) denote re- 
spectively the pseudo-remainder and the pseudo-quotient of F w.r.t. G in Xk- For any 
triangular set T = [T\, . . . ,T r ] and polynomial set V in lC[x], define prem(F, T) := 
prem(- • • prem(F, T r ,lv(T r )), . . . , Ti, lv(Ti)) and premCP, T) ■= {prem(P, T) : P E V}. 

Definition 8. A (weak) ascending set C in IC[x] is called a (weak) characteristic set of 
a non-empty polynomial set V C K,[x] if C C (V) and prem^C) = {0}. 

For characteristic sets, the following property is fundamental. 



Proposition 9 (|Wulll986l) . LetC = [Ci, . . . , C r ] be any characteristic set of a polynomial 
set V C ]C[x] and 

Ii=im(Ci), V i = PU{I i }, X={h,...,I r }. 

Then 

Zero(C/Z) C Zcro(^) C Zero(C), 

(1) 



Zero(T') = Zero(C/J) U |J Zero(Pi). 



3. Algorithmic Scheme for Computing Generalized Characteristic Sets 

In this section we generalize the concept of characteristic set and give a proposition 
to show that the generalized concept preserves the basic property (Proposition 9) of the 
original one. For computing generalized characteristic sets, we propose a new algorithmic 
scheme, in which the set of reductions, the termination condition, and strategies for find- 
ing reduction polynomials can be viewed as placeholders. The scheme will be instantiated 
as concrete algorithms in Section 5. 

Definition 10. For any non-empty polynomial set V in JC[x], an ascending set C C K.[x] 
is called a generalized characteristic set of V if 

(a) CC(P>; 

(b) there exists a polynomial set Q C K,\x] such that (Q) = (V) and prem(Q, C) = {0}. 

It is easy to see that characteristic sets are indeed a special kind of generalized char- 
acteristic sets. The zero relations in (1) also hold for generalized characteristic sets. 

Proposition 11. Let C = [Ci, . . . , C r ] be any generalized characteristic set of a polyno- 
mial set V C K\x\ and 

Ii=im(Ci), P i = PU{I i }, X={h,...,I r }. 

Then the zero relations in (1) hold. 

Proof. By Definition 10, there exists a Q C )C[x] such that (V) = (Q) and prem(Q,C) = 
{0}. It is easy to know that C is a characteristic set of Q. Bearing in mind Proposition 9 
and Zero('P) = Zero(Q), one sees clearly the truth of the proposition. □ 
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Input: T — non-empty polynomial set in /C[a;]; 

D — set of admissible reductions in K\x\. 
Output: A — generalized characteristic set of T . 

o.i Q := 1Z := J-: 

0.2 while 72. 7^ do 

o.3 [A,B] := MedSet(£,D); 

0.4 if A is a contradictory ascending set then 
0.5 K:= 0; 

o.6 else 

0.7 1Z:=prem(B\A,A)\{0}; 
o.s g --guAUK; 

0.9 end 

o.io end 



Now we can describe the main algorithm NewCharSet (Algorithm 1), which computes 
a generalized characteristic set of any given non-empty polynomial set in IC[x]. 



Algorithm 1 is similar to the algorithm GenCharSet presented in IWanj ( 2001bl ). The 



difference is that MedSet here outputs not only a medial set A of the polynomial set 
g, but also another basis B of the ideal (Q). Properly storing in B the information 
from computing A, one may check more efficiently whether A is a characteristic set, 
which usually takes a large proportion of the total computing time. Before describing 
the concrete steps of MedSet, we need some additional definitions such as admissible 
reduction. 

Definition 12. Let P,Q E K[x] \ {0} and <i ox be the lexicographic order in K.[x] such 
that x\ < ■ ■ ■ < x n - Denote P < Q or Q > P if 

(a) ht(P) <i cx ht(Q), or 

(b) ht(P) = ht(Q) and P - hm(P) <Q- hm(Q), 

where ht(P) and hm(P) stand for the heading term and the heading monomial of P 
under <i cx respectively. Set < P for any P e JC{x] \ {0}. 

If neither P < Q nor P > Q, we write P w Q. The notation P ^ Q means that P > Q 
or P ~ Q, and P ^ Q is similarly defined. 

The partial order < is a refinement of -< (Definition 3) . More precisely, for any P,Q € 
JC[x\ \ {0}, P -< Q implies that P < Q. On the contrary, P < Q implies only P ^ Q, and 
P -< Q does not necessarily hold. Consider for example the polynomials P — y 3 + x and 
Q = y 3 + x 2 in fC[x, y]. By definition, P < Q and P ~ Q. 

Lemma 13. For any polynomial sequence Pl ~ P2 « P3 ~ ■ • ■ ? there exists an integer 
m such that for any i > m, Pi w P m . 
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Proof. From Pi P 2 P3 <J • • • , one can obtain the seq uence ht(Pi) > ht(P2) > 
ht(Ps) > • ■ • . In view of the property of term orders (see, e.g.. lCox et al. 19971 Lemma 2 
in Chapter 2), there exists an integer s such that for any i > s, ht(Pi) = ht(P s ). 

Let P/ = Pi — hm(P,). We have a new sequence P' s P s ' +1 <; P s ' +2 ^ • • • . Similarly, 
there exists an integer t (> s) such that for any i > t, ht(P/) = ht(P/). The rest can be 
done in the same manner. As the number of terms of any polynomial is limited, we will 
finally obtain an integer m such that for any i > m, all P; have the same set of terms, 
which means Pi ~ P m . □ 



Let 25 be an operation with two polynomials P,Q G lC[x]\K, as its input and an ordered 
set of two polynomials in IC[x] as its output, and the latter is denoted by Rem(P, Q, 25). 
The definition of admissible reduction, which is somewhat abstract, is given below. In the 
next section, we will discuss a number of concrete reductions used in our experiments. 

Definition 14. Let [Ri,R 2 ] = Rem(P,Q,25) for any P, Q G K[x] \ K. The operation 25 
is called an admissible reduction (or reduction for short if there is no confusion) in K,[x] 
if Pi,P 2 G (P,Q). 

Suppose that 25 is an admissible reduction in K\x\. We say that P is Id-reducible w.r.t. 
Q if P < Pi and Q ;$ R2; otherwise, P is D-reduced w.r.t. Q. Moreover, P and Q are 
called the reductend and the reductor respectively, and Rem(P, Q, 25) the 25 reduction-rest 
of P by Q. 

The subalgorithm MedSet is given as Algorithm 2, where cond can be replaced by any 
Boolean expression that guarantees the termination of the while loop. 

The operation Find3R(«4, D) chooses polynomials P (reductend) and Q (reductor) from 
A and 25 (reduction) from KS such that P is 25-reducible w.r.t. Q. If Find3R(^l, D) finds 
appropriate P, Q,25, then it returns this triple; otherwise, it returns []. 

The function RemCh is an extension of Rem. RemCh(P, Q, 25) returns not only the 25 
reduction-rest TZ of P by Q, but also a Boolean value b. When b is true, the reduction-rest 
TZ must satisfy the condition P, Q G (TZ) . The computation of b depends on the concrete 
admissible reduction 25, which will be discussed in the next section. The if block (lines 
0.10-0.12) is designed to store information acquired in the reduction process, which may 
be used by NewCharSet to check more efficiently whether the output Ai is a characteristic 
set. 

In algorithm MedSet, one may view cond, 25 and Find3R as placeholders, which will 
be instantiated in Section 5. The correctness and termination of Algorithms 2 and 1 are 
proved as follows. 



Proof. (Algorithm 2) Correctness. For the statement [TZ,b] := RemCh(P, Q, 25) in line 
0.4, we have TZ C (P, Q) by Definition 14. It is thus easy to verify that A C (J 7 ) always 
holds during the running of the algorithm, which means M. C (J 7 ) . Obviously, the ranking 
of M. is not higher than the ranking of any basic set of J 7 ; hence is a medial set of T 
by definition. 

Let Qi be the initial value of Q in the ith while loop. Obviously, (Q\) = (J 7 ). Suppose 
that (Qi) = (J-), and we assert that (Gi+i) = (J 7 ) as follows. 

Consider the statement [TZ, b] := RemCh(P, Q, 25) in line 0.4 of the ith while loop. If 
the Boolean expression (P, Q G Q and b) equals false, then Gi+i = Qi- Hence (Qi+\) = 
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Input: T — non-empty polynomial set in /C[a;]; 

D — set of admissible reductions in K\x]. 
Output: M. — medial set of J 7 ; 

Q — polynomial set satisfying (Q) = {J 7 }. 

o.i A := Q := J-: 
0.2 while cond do 

0.3 [P,Q,S] := Find3R(^,B); 

o.4 [K,b] := RemCh(P,Q,£>); 

0.5 if 1Z contains a non-zero constant then 

o.6 A := {1}; Q := {1}; break; 

0.7 else 

o.s A:=A\{P,Q}UK\{Q}; 
0.9 end 

o.io if P, Q G Q and b then 

G :=g\{P,Q}UK; 

0.12 end 
0.13 end 

o.i4 M := BasSet(.4U J 7 ); 



(J 7 ). Now assume that (P, Q G Q and b) equals true. Then line 0.11 is executed, and we 
have Qi+i = Qi \ {P, Q} U 1Z. As D is an admissible reduction, we have 1Z C (P, Q). By 
the assumption on b returned by RemCh, one has P,Q& (1Z). Thus (P, Q) = (1Z), which 
implies that (Q i+1 ) = (Qi \ {P,Q} U 1Z). From the above, {Q i+1 ) = (J"), so (Q) = (J 7 ) 
always holds. 

Termination. It is obvious by the assumption on cond. □ 

Proof. (Algorithm 1) Correctness. From the properties of MedSet's output, we know 
that (B) = (Q) in line 0.3, and A is a medial set of Q. It follows that A C (Q). According 
to the pseudo-division formula, for any polynomial R in 7Z of line 0.7, there exist B S S\^4 
and Gi G K.[x] such that 

(Q ini(^) 8 *) • B = GiAi + R, 

i i 

where each s, is a non-negative integer and Aj G A. Hence 1Z C (Q) always holds during 
the running of the algorithm. Thus the ideal generated by Q remains the same after the 
assignment in line 0.8, which means that (Q) = (T) always holds. 

On the other hand, piem(Q,A) = {0} when the while loop terminates. Hence A is a 
generalized characteristic set of J- by Definition 10. 

Termination. We use A% and Qi to denote the values of A and Q respectively in the 
ith while loop after executing line 0.3. Recalling the properties of MedSet, we know that 
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each Ai is a medial set of Qi. By Proposition 7, one can obtain the sequence A\ >~ A2 >- 
A3 >- ■ ■ • of triangular sets, which should be finite by Proposition 5. Thus the algorithm 
terminates. □ 

4. Concrete Admissible Reductions 

In this section we introduce and discuss several concrete admissible reductions. 
®UG {univariate GCD reduction): define Rem(P, Q, Sue) : — 

{[0, gcd(P, Q, x q )], if P, Q are univariate polynomials in x q ; 
[P,Q], otherwise. 

In the above definition, x q is some variable in x. It is easy to verify that S?tjg is an 
admissible reduction, and P is Sue-reducible w.r.t. Q if and only if P,Q are univariate 
polynomials in the same variable. 

Pseudo-division, used frequently in triangular decomposition, is also an admissible 
reduction. 

Dp (pseudo-division reduction): define Rem(P, Q, Dp) := [prcm(P, Q, lv(Q)), Q]. We 
have the following pseudo-division formula 

ini(Q) 6 P = pquo(P, Q, lv(Q)) Q + prem(P, Q, lv(Q)), (2) 

where s is a non-negative integer. Thus Rem(P, Q, Dp) C (P,Q), which means that Dp 
is an admissible reduction in K\x\. 

Proposition 15. P is Dp-reducible w.r.t. Q if and only if P is reducible w.r.t. Q. 

Proof. Let x q = fv(Q). 

(=>) Suppose that P is reduced w.r.t. Q. According to the definition of pseudo-division, 
we have prem(P, Q, x q ) = P. Thus P is 2)p-reduced w.r.t. Q. 

(<=) Suppose that P is reducible w.r.t. Q. Then either 

prem(P, Q, x q ) = or deg(prem(P, Q, x g ), x q ) < ldeg(Q). 

If prem(P, Q, x q ) = 0, then P is obviously £>p-rcducible w.r.t. Q. Now assume that 
prem(P, Q, x q ) ^ 0. By analyzing the algorithm of pseudo-division (see, e.g. jMishralll993l 
PseudoDivisionRec), one may find that for any i(q < i <n), deg(prem(P, Q, x q ),x t ) < 
deg(P, Xi). As P is reducible w.r.t. Q, 

deg(prem(P,Q,a; g ),x ( j) < ldeg(Q) < deg(P,x g ). 

Hence prem(P, Q, x q ) < P. It follows that P is S)p-reducible w.r.t. Q. □ 

Pseudo-division can be done step by step using the following operation. 

Definition 16. Let P, Q e JC[x] \ JC with x q = lv(Q), I = ini(Q), and J = \c(P,x q ) 
and suppose that P is reducible w.r.t. Q. Then we can perform the following one-step 
pseudo-division: 

R :=FP- GQ x ^s(P,x q )-\d cg (Q) ^ ^ 

where F = lcm(/, J)/ J and G = lcm(/, J)/ 1. R is called the one-step pseudo-remainder 
of P w.r.t. Q and denoted by stprem(P, Q). 
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One can see that pseudo-division is recursive application of one-step pseudo-divisions 
and thus may immediately lead to big superfluous factors of the pseudo-remainder. In 
contrast, it is easier to control the selection and size of polynomials for one-step pseudo- 
division, which may result in smaller reduction-rests when combined with other reduc- 
tions. 

®sp {one-step pseudo- division reduction): define Rem(P, Q 1 Dsp) := 

{[stprem(P, Q), Q], if P is reducible w.r.t. Q; 
[P, Q] , otherwise. 

If P is reducible w.r.t. Q, then stprem(P, Q) € (P, Q) by (3), and hence Rem(P, Q, Dsp) C 
(P, Q); otherwise, Rem(P, Q, 25sp) = [P, Q] C (P, Q) is obvious. Therefore £»sp is also an 
admissible reduction by definition. 

Proposition 17. P is Ssp -reducible w.r.t. Q if and only if P is reducible w.r.t. Q. 

Proof. If P is reduced w.r.t. Q, then Rem(P, Q, £>sp) = [PQL an d hence P is 
£>sp-rcduced w.r.t. Q. 

(<=) Suppose that P is reducible w.r.t. Q. Then Rem(P, Q,£sp) = [stprem(P, Q), Qj. 
If stprem(P, Q) = 0, then P is obviously 2)sp-rcducible w.r.t. Q. Now assume that 
stprcm(P, Q) ^ and consider (3). We know that 

F = lcm(/, J) /J = 1/ gcd(J, J) G fc[a?i, . . . , 

Thus for any i (q < i <n), deg(FP,Xi) = deg(P, Xi). Furthermore, 

deg(G0^ cg ( p '^)- ldcs ( Q ),x J ) = deg(G,Xi) = dcg( J/ gcd(7, J), x r ) < deg(P,x I ). 

Hence deg(stprem(P, Q), < deg(P,Xi). Moreover, we have 

deg(FP lXq ) = deg(GQxf g(p ' :!: ' ) - ldos(Q) ,x (3 ), 

which implies that deg(stprem(P, Q),x q ) < deg(P, x q ). It follows that stprem(P, Q) < P, 
so P is £>sp-reducible w.r.t. Q. □ 

The division operation, which is used in the computation of Grobner bases, can also be 
viewed as an admissible reduction. Instead of introducing the division reduction directly, 
we discuss the one-step division reduction first and then use it to recursively define the 
former. 

2)sd (one-step division reduction): define Rem(P, Q, S)sd) •= 

[P — Mj ht(Q) • Q, Q), if there exists a monomial M of P such that ht(<5) | M; 
[P, Q], otherwise. 

It is easy to verify that 2Dsd is an admissible reduction, and P is 2DsD-reduciblc w.r.t. Q 
if and only if there exists a monomial of P which can be divided by ht(Q). 

Bearing in mind the relation between 2Dsp and Dp, we can define the division re- 
duction £>d as follows. First, set Rq := P and Sq := Q. Then recursively compute 
[Pi+i, := Rem(Pj, 5^, Ssd), until an integer m is found such that R m +i = R m - 

Define Rem(P, Q, Sp) :~ [R m , S m ], where S m = Q. Similar to S3sdj S)d is an admissible 
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reduction, and P is 2) d -reducible w.r.t. Q if and only if there exists a monomial of P 
which can be divided by ht(Q). 

By computing subresultant PRS, one can also design a useful admissible reduction. 
We gi ve the definition here and the reader may refer to Ivon zur Gathen and Luckind 
(|2003l ) for more details. Suppose that lv(P) = lv(Q), ldeg(P) > ldeg(Q) and treat P,Q 
as univariate polynomials in x q = lv(Q) with coefficients in K[x\, . . . , x q -{\. Define the 
subresultant PRS of P and Q w.r.t. x q to be 



Pi '■— Pi Pi '■— Q, Pi+2 '■— 



prem(P i ,P i+ i,a; ( j 



ii+2 



1 < i < r - 2, 



where 



Qi 

H, 



-1, 



deg(Pi,x,j), i=l,...,r, 
-\c(P- 2 ,x q )H?*- 2 - d ' 

(-lc(P_ 2 ,X g )) d - 3 -*- 2 



i=4,..., 



"gJJ i—1 ' — ' ' ' ' ' ' 

and prem(P r _i, P r , x q ) = 0. The well-known Subresultant Chain Theorem (Mishra 19931 
Theorem 7.9.1) indicates the relation between the subresultant PRS and the Euclidean 
PRS of P and Q: Pi is similar to the element of the same degree in the Euclidean PRS. 
Furthermore, the former may have smaller superfluous factors. 

We use res(P, Q, x q ) to denote the resultant of P and Q w.r.t. x a . If res(P, Q, x q ) = 0, 
then P r is a greatest common divisor of P and Q w.r.t. x q (|Mishrall 19931 Corollary 7.7.9). 
If res(P, Q, x q ) ^ 0, then cls(P r ) < cls(P) and cls(P r _i) = cls(P). Thus P r _i is a greatest 
common divisor of P and Q w.r.t. x q under the condition P r = ( PI ; on the other hand, 
under the condition P r ^ 0, P r would be a greatest common divisor of P and Q. Hence 
we can introduce the following reduction. 

J)sc {subresultant PRS reduction): define Rem(P, Q, Ssc) := 

[0, P r ], if lv(P) = lv(Q), ldeg(P) > ldeg(Q) and rcs(P, Q, lv(Q)) = 0; 
[P r ,p._i], if lv(P) = lv(Q),ldeg(P) > ldeg(Q) and res(P, Q, lv(Q)) ^ 0; 
[P, Q] , otherwise. 



Example 4.1. Let P 

of P and Q in x is a x 2 + 



£> sc , Rem(P, Q,D S c) = {d 2 c- edb + ae 2 ,dx + e]. 



ax 2 + bx + c and Q = dx + e. Then the subresultant PRS 
b x + c, dx + e, <i 2 c — e<i& + ae 2 . According to the definition of 



According to Mishra! (|l993i Lemma 7.7.4), there exist Aj, Bj £ K[x] such that AjP + 
BjQ = Pj, i.e., Pj G (P, Q). Thus Dsc is an admissible reduction. In addition, it is easy 
to verify that P is 2)sc- rc duciblc w.r.t. Q if and only if lv(P) = lv(Q) and ldcg(P) > 
ldeg(Q). If ldcg(P) < ldeg(Q), then consider Rem(Q, P, D sc ). 

The standard reduction used in Ritt-Wu's algorithm CharSet and its variants is pseudo- 
division only, which may lead to superfluous factors and the swell of polynomial coef- 
ficients. This is one of the main causes that degrades the performance of Ritt-Wu's 



An accurate description requires the use of evaluation homomorphism, see Mishra ( 1993|, Section 7.3.1) 
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algorithm in many cases. In the algorithm NewCharSet (line 0.7) other kinds of admissi- 
ble reductions may also be used. The condition of admissible reductions is quite weak, 
which makes our algorithm clearly more flexible than other existing ones for computing 
characteristic sets. 

The remaining problem is the computation of b in the step [lZ,b] := RemCh(P, Q, 33) 
of Algorithm 2, where P is ©-reducible w.r.t. Q. 

Let [i?i,i?2] = P- For 33tjg, we have Pi = and R 2 — gcd(P,Q,x q ), which implies 
that (P, Q) C (Pi,P 2 ). For 33 S d and 33 D , R 2 — Q and the formula R l = P - AQ always 
holds, where A is a polynomial in JC[x]. Thus it is obvious that P £ (Ri,Q) = (Ri,R 2 ) 
and Q G (Pi, R 2 ). Hence for any of 33tjg, ©sd, and 33d, we can let RemCh return true 
as the value of b. 

For 33p, recall the pseudo-division formula (2). Provided that mi(Q) s is a non-zero 
constant, we also have P £ (Ri,Q) = (Ri,R 2 ). Thus for the pseudo-division reduction, 
we set b to be the value of the Boolean expression (ini(Q) £ K or s = 0). Similarly, 
consider (3) for 33sp and set b to be the value of the Boolean expression (F £ JC). 

Finally for 33sc, s °t b = false permanently. 

5. A Sample Subalgorithm 

In this section we provide details about the algorithm MedSet. The assignment of the 
Boolean expression cond is first considered and the termination of Algorithm 2 guaranteed 
by cond is proved. Then the basic idea and a sample algorithm for the operation Find3R 
are presented. We conclude this section with discussions about other possible strategies 
for Find3R. 

There are several ways to assign cond. For example, one can set cond = false. In 
this case, the while loop in MedSet never starts and NewCharSet is identical to Ritt-Wu's 
algorithm. Another example is to set a counter i of the while loop and assign an inequality 
to cond, say i < 50; then the while loop will run 50 times. The termination of MedSet in 
both cases is obvious. In our implementation, we let cond be Find3R(^4, D) ^ [], which 
means that the while loop repeats until there is no triple [P, Q, 33] (P, Q £ A, P ^ Q, 33 £ 
D) such that P is 33-rcduciblc w.r.t. Q. For this case, the termination of MedSet is proved 
as follows. 

For any polynomial Pi in J 7 , if we perform [R 2 ,Q 2 ] := Rem(Pi, Qi, 33), where Qi 
is chosen by Find3R such that Pi is 33-reducible w.r.t. Qi, then R 2 can be viewed as 
the successor of Pi with Pi > P 2 ; if we perform [P 2 ,R 2 ] := Rem(Pi, Pi, 33'), where Pi 
is chosen by Find3R such that Pi is 33'-reduciblc w.r.t. Pi, then R 2 can be viewed as 
the successor of Pi with Pi <; R 2 . As the algorithm runs, the successor P 3 of R 2 and 
then the successor P4 of P3 may be produced in the same way. Thus one may obtain a 
polynomial sequence Pi ^ P2 ~ P3 ~ ■ • • . By Lemma 13, there exists an integer m such 
that P m « Pm+i ~ • • ■ • As this is the case for every polynomial in T , the operation 
Find3R will find no triple in the end. Hence the algorithm terminates. 

Now we explain the basic idea of the algorithm Find3R. The set of selected admissible 
reductions is D = {33tth. 33gp, 33sc, 33qp}l 2 I. There may exist several triples and/or sev- 
eral admissible reductions between a pair of polynomials. We would like to select better 



2 £>d an< i 35p are omitted because they can be treated as recursion of £>sd and X>sp respectively. 
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triples first. Therefore, a way to measure the goodness of triples should be defined. In 
addition, an order of the admissible reductions in D is also needed. We give a sample 
order on triples as follows. 

First of all, it is natural that a triple with S)ug is better than other triples without it. 
In fact, several univariate pseudo-remainders may be produced in line 1.7 of the algorithm 
NewCharSet. It is easier to deal with the Dug reduction-rest of such pseudo-remainders 
in the next while loop than to deal with them immediately Next, a triple with T>sd 
should be better than a triple with 2)sp or Ssc- There are two reasons for this. 

(1) For two triples [P,Q,Dsu] and [P,Q,D S p], P is both S SD -reducible and £>sp- 
reducible w.r.t. Q. It is easy to see that computing the ®sd reduction-rest R of 
P by Q requires less multiplications than computing the Ssp reduction-rest R' 
and in general R is of smaller size than R' . Therefore, [P,Q,®sd] is better than 

[P,Q,Qsp]. 

(2) For two triples [P, Q, D SD ] and [P', Q', D S p] with P ^ P' or Q ^ Q' , the former is 
better because P is in (Q,R), while in general P' is not in (Q',R'), where R and 
R' are the 2)sd reduction-rest of P by Q and the £>sp reduction-rest of P' by Q' 
respectively. 

Then, [P, Q, £>sc] may be better than [P, Q, £>sp] because the reduction-rest of the former 
generally involves coefficients of smaller size than that of the latter. Finally, a triple 
[P, Q, S] might be better than [P\ Q', D] (2 e D) if P > P' or Q < Q' . This observation 
is based on the way of reduction in Ritt-Wu's algorithm: the polynomial of maximal class 
and maximal degree is first reduced by a polynomial of minimal degree. To summarize, 
we provide the formal definition of an order on the triples: [P, Q, 5)] < [P' ,Q' ,&] if 
D < ID', or D = £>' and P < P', or D = £>' and P » P' and Q > Q'. The admissible 
reductions in D are ordered as £>ug > ®sd > ®sc > ®sp- 

Taking as input a polynomial set J 7 and a set D of selected admissible reductions under 
the above order, Find3R will do the following steps to select a better triple or output if 
there is none. 

1. Check if there exist triples [P, Q,33tjg] (P,Q G J 7 , P ^ Q). If there is only one, 
then return it; if there are many, then return the one in which P has maximal class 
and maximal leading degree and Q has fewest terms and minimal leading degree. 

2. Sort the polynomials in T increasingly w.r.t. the partial order <. Start with a 
polynomial P of highest order and check if there exists any polynomial Q such that 
P is SsD-rcducible w.r.t. Q. If there is only one Q, then return [P, Q, £>sd] ; if there 
are many, then select Q from 

{F e T : P is 2)sD-reducible w.r.t. P} 

which has fewest terms and minimal leading degree and return [P, Q, Qsd]] other- 
wise, consider another P of lower order and check again. If there is no such triple, 
then go to the next step. 

3. Start with a polynomial P of highest order again and check if there exists any Q 
such that P is Ssc-rcducible w.r.t. Q. The process is similar to step 2. 

4. Start with a polynomial P of highest order again and check if there exists any Q 
such that P is Ssp-reducible w.r.t. Q. The process is similar to step 2. 

5. If there is no triple at all, then output 0. 
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Input: JF — set of polynomials in K\x\\ 

D= [D ug ,S)sd,S)sc,25sp]- 
Output: or [P,Q,D] such that P is ©-reducible w.r.t. Q. 

o.i S := 0; 

0.2 for i = n to lwhile |5| < 2 do 

o.3 S := {P G T : P G K,[xi]}; 
0.4 end 

0.5 if 3 P, Q G S, P 7^ Q, such that P is Qua-reducible w.r.t. Q then 

return [P, Q, Sue] such that P e S has maximal degree, Q G S \ {P} has 
fewest terms and minimal degree, and P is 2)uG _re ducible w.r.t. Q; 

end 

[Pi, . . . ,P r ] := J- (with P; sorted increasingly w.r.t. <); 
for i — r to 2 do 

Q := {Q G P\ {P} : P is S SD -reducible w.r.t. Q}; 
if Q ^ then 

choose Q G Q with fewest terms and minimal leading degree; 
return [P,Q,D S d]; 

end 

0.15 end 

o.i6 for i = r to 2 do 

Q := {Q G J"\{Pi} : P is D SC -reduciblc w.r.t. Q}; 
if Q 7^ then 

choose Q G Q with fewest terms and minimal leading degree; 
return [P,Q,S SC ]; 

end 

end 

for i — r to 2 do 

Q := {Q G J"\{Pi} : P is S S p-reduciblc w.r.t. Q}; 
if Q / then 

choose Q G Q with fewest terms and minimal leading degree; 
return [P,Q,£>sp]; 

end 

end 
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The above process is described formally as Algorithm 3. Its correctness follows from 
the above analysis and its termination is obvious. If "P* is £>sp-reducible w.r.t. Q" in line 
3.24 is replaced by "ini(Pj) is Ssp-reduciblc w.r.t. Q" , then MedSet computes a weak- 
ascending set and therefore NewCharSet outputs a weak-characteristic set. 

Note that the design of Find3R is flexible and can be made more technical and compre- 
hensive. There are several ways to improve Find3R. One may introduce other admissible 
reductions by computing Bczout resultants or Grobner bases, or using other techniques 
(e.g., from linear algebra). If the input set of Find3R contains polynomials of special form 
or structure, then new reduction strategies may be adopted. For example, if there is one 
univariate polynomial of low degree or few (say one or two) terms, then one should use 
this polynomial first to reduce (and simplify) other polynomials. The ordering used in 
line 3.8 is crucial for the efficiency of MedSet. One may sort polynomials w.r.t. different 
orderings depending on the admissible reductions. 



6. Example and Preliminary Experiments 

In this section we present first an illustrative example to compare the outputs of 
different algorithms and to show how the algorithm NewCharSet works and then our 
experimental results on the performance of NewCharSet. In what follows, an index tuple 

[[deg(P,xi),.. .,deg(P,x n )],nops(P),hm(P),m] 

is used to characterize a polynomial P £ K.[x] \ fC, where nops(P) denotes the number of 
terms of (expanded) P, hm(P) the heading monomial of P, and m the maximal number 
of digits of the integer coefficients of P. 



Example 6.1 (|Wanj20o4 Epsilon-A14). Let T= {Fx,F 2 ,F 3 } C Q[w,x,y,z], where 



Pi = x 2 + y 2 + z 2 — w 2 , P2 = xy + z 2 — 1, P3 = xyz — x 2 — y 2 — z + 1, 

with w < x < y < z. The output of CharSelT^I consists of three polynomials with index 
tuples 

[[520, 42, 0, 0], 5450, w 508 a; 42 , 212] , 

[[6,30,l,0],93,x 29 y,3], (4) 

[[520, 40, 0, 1], 10390, w 510 x 40 z, 212] , 

while the output of NewCharSet also consists of three polynomials yet with index tuples 

[[8, 12, 0,0], 23, z 12 , 2], [[4,6,l,0],12,™ 2 x 3 y,l], [[4, 6, 0, 1], 17, x 6 z, l] . (5) 

The Grobner basis of T w.r.t. the lexicographic order (determined by w < x < y < z) 
contains 9 polynomials, of which one of the biggest has index tuple [[12, 11, 1, 1], 41, yz, 5] . 

Now let us show how NewCharSet works. For brevity, the admissible reductions Dsn 
and £>sp are renamed £>d and Dp respectively. After the initialization of Q and A in line 
1.1 of NewCharSet, the first while loop (line 1.2) starts and MedSet(P, B) is called. Line 
2.1 initializes the values of A and Q. Then it is checked whether the Boolean expression 



Corresponding to the Maple command charse^J 7 , [w, x, y, z]) with the Epsilon package I Wana 2004 ) 



15 



Find3R(.A, D) ^ [] is true or false. One sees clearly that F3 < F2 < F\ and there are 
three triples 

[F 2 ,F 3 ,!Dd], [Fi.Fg.^D], [F x ,F 2 ,®v\. 
Therefore the first while loop (line 2.2) starts. Find3R(^4., B) outputs the third triple 
and RemCh(Fi, F2, £>d) then returns [Rem (F\, F2, 23d), true]. Since the reduction-rest 
contains no constant, A and Q are updated and the second while loop (line 2.2) starts. 
After the second while loop, Q is updated to be 

{y 2 — xy + x 2 — w 2 + 1, xyz — z — xy — w 2 + 2, z 2 + xy — 1} (6) 

and remains unchanged afterwards. After the 7th while loop, A is updated to be an 
ascending set with index tuples of its polynomials shown in (5) and Find3R(^4, D) = []. 
So MedSet returns a basic set of A U T (which is A) as well as (6). It can be verified by 
line 1.7 that the pseudo-remainders of the polynomials in (6) w.r.t. A are all zero. Note 
that it is easier to compute prem(C7 \ A, A) than prem(J-"\ A, A) because one polynomial 
in (6) is of smaller class than the polynomials in T . Therefore, A with the index tuples 
in (5) for its three polynomials is a characteristic set of T . 

We have implemented the algorithm NewCharSet described in Section 3 using the 
Epsilon library. There are two versions of it: NewCharSetw outputs weak-characteristic 
sets and NewCharSet outputs characteristic sets. We have compared the performances 
of NewCharSet, NewCharSetw, CharSet, and CharSetwFTas well as the Grobner basis al- 
gorithm Groebneff]- The timings (in CPU seconds) given in Tabic 1 arc for our imple- 
mentation running in Maple 14 under Windows 7 on a laptop Intel Core 2 Duo T6670 
Processor 2.20 GHz with 3 GB of memory. The asterisk * indicates that the computation 
is out of memory. Table 2 shows the size of the test examples (including the number of 
variables and the total degree of each polynomial) and the number of polynomials in the 
output of each algorithm. The entry [3, 3, 4] in the table means that example 1 consists of 
three polynomials whose total degree are 3, 3, 4 respectively. The test examples are taken 



from|http: / /www. symbolicdata.org/ and http://www-salsa.lip6.fr/" wang/cpsilon/. It is 



easy to see that NewCharSet is more efficient than CharSet and CharSetw for most of the 
examples and is comparable with Groebner. 

Now let us comment on the empirical results shown in Table 1. For examples 3, 5, 
8, 10, 15, and 16, NewCharSet is slower than CharSet, and for examples 9, 13, 17, and 
18, neither of them can get any result within 1000 seconds. For the 8 other examples 
(notably, examples 4 and 14), NewCharSet is always faster than CharSet. Among the first 
four algorithms, NewCharSetw is the most efficient (for all but one test example, i.e., 
example 15). There are five examples for which CharSet cannot get any result within 
1000 seconds. The Groebner algorithm is slower than other algorithms for examples 1, 3, 
7, 14, and 15, but considerably faster for examples 9 and 18. There are two examples for 
which Groebner cannot complete the computation due to the lack of memory. 

There are several reasons for NewCharSet to be more efficient than CharSet. First of all, 
it is the polynomials in the output Q of MedSet (instead of J 7 ) whose pseudo-remainders 



4 By CharSet and Ch a rSetw we me an the algorithms implemented in Epsilon which adopted optional 
strategies from Wane 1 19921 . 2001a ) for speeding up the computation of characteristic sets rather than 



Ritt-Wu's original algorithm. 

Corresponding to the Maple command Groebner [Basis] (J 7 , plex(x„, . . . ,xi)). 
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Table 1. Timings (in seconds) for five algorithms 



No 


Name 


CharSetw 


CharSet 


NewCharSetw 


NewCharSet 


Groebner 


1 


DiscrC2 


0.016 


0.016 


0.016 


0.016 


>1000 


2 


Epsilon-A14 


>1000 


16.302 


0.016 


0.031 


0.156 


3 


Choul56-l 


0.218 


0.093 


0.047 


0.156 


19.812 


4 


Trinksl 


>1000 


>1000 


0.063 


0.063 


0.047 


5 


ZeroDiml4 


34.398 


0.795 


0.156 


7.909 


0.172 


6 


Epsilon-A16 


761.145 


8.128 


0.390 


3.073 


0.266 


7 


Schiclcl 


1.576 


1.482 


0.406 


0.375 


5.491 


8 


Epsilon-A25 


0.172 


0.078 


0.546 


7.301 


1.029 


9 


Cyclic5 


>1000 


>1000 


0.577 


>1000 


0.062 


10 


Feci 


o7o.0ol 


2.777 


1.326 


5.273 


0.156 


11 


Weispfenning94 


>1000 


33.680 


2.340 


7.238 


0.921 


12 


Stcidcl2 


>1000 


49.312 


5.475 


15.678 


0.905 


13 


Epsilon-A30 


>1000 


>1000 


7.941 


>1000 


2.450 


1 1 


Fcitemcin 


>1000 


>1000 


9.797 


24.429 


411.546 


15 


Sym3-5 


23.868 


1.623 


10.733 


10.483 


35.865 


16 


Epsilon-A3 


160.150 


67.174 


17.035 


261.192 


* 


17 


Wu90 


>1000 


>1000 


44.835 


>1000 


* 


18 


Cyclic6 


>1000 


>1000 


>1000 


>1000 


0.374 



w.r.t. the output (weak-) medial set M. of MedSet are computed. This may reduce the 
cost of zero pseudo-remainder verification in most cases because the polynomials in Q are 
"closer" (than those in F) to the polynomials in M.. For example, if Q in MedSet is fixed to 
be J 7 , then NewCharSet takes 42.307 seconds (vs 15.678 seconds when Q remains updated) 
for example 12 and gets no result within 500 seconds (vs 24.429 seconds when Q remains 
updated) for example 14. However, keeping Q updated may also lower the efficiency of 
the algorithm for some examples (including three of our test examples) . How to choose 
Q to improve the efficiency is a question that remains for further investigation. 

The use of admissible reductions allows us to effectively control the swell of poly- 
nomial coefficients as well as degrees. As we have seen, the maximal numbers of digits 
of coefficients in (4) are much bigger than those in (5). More comparisons are given in 
Table 6 in the appendix. The ordering < used in line 3.8 can also be replaced by other 
orderings. According to our experiments, there may be fewer while loops in NewCharSet 
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Table 2. Numbers of polynomials in the outputs of five algorithms 



No No of vars Total degree list CharSetw CharSet NewCharSetw NewCharSet Groebner 



1 


12 


[3, 3, 4] 


3 


3 


3 


3 




2 


4 


[1, 2, 3, 4] 




3 


3 


3 


9 


3 


1 


[2, 2, 2, 2] 


4 


i 


4 


4 


49 


4 


6 


[1, 1, 2, 2, 2, 3] 






6 


6 


6 


5 


4 


[1, 2, 3, 4] 


1 


1 


4 


1 


1 


6 


4 


[4, 3, 2, 4] 


4 


i 


1 


4 


1 


7 


2 


[6, 11] 


2 


2 


2 


2 


3 


8 


8 


[2, 2, 2, 2] 


1 


1 


4 


4 


21 


9 


5 


[1, 2, 3, 4, 5] 






5 




11 


10 


4 


[2, 3, 4, 4] 


•1 


i 


1 


4 


5 


11 


3 


[4, 5, 5] 




3 


3 


3 


3 


12 


3 


[4, 5, 5] 




3 


3 


3 


3 


13 


10 


[3, 3, 3, 3, 3, 3, 3, 3, 3, 3] 






10 




10 


14 


1 


[3, 5, 5] 






3 


3 


108 


15 


3 


[6, 6, 6] 


3 


3 


3 


3 


3 


16 


4 


[4, 4, 4, 5] 


4 


i 


4 


4 




17 


4 


[3, 3, 3, 4] 






4 






18 


6 


[1, 2, 3, 4, 5, 6] 
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if polynomials are sorted w.r.t. the lexicographic order of their degree tuples, i.e., P is 
lower than Q if 

[deg(P,xi), . . . , deg(P, x n )] < lox [deg(Q, xi), . . . , deg(Q, x„)]. 

For example, let P = yz + x 2 z + 1 and Q = xz + z + y 2 with x < y < z; then P > Q 
according to Definition 12. However, the above ordering indicates that P is lower than 
Q. One can see that the degree tuples of P and Q are [2, 1, 1] and [1, 2, 1] resp ectively, so 



[2, 1 11 <i cx [1, 2, 1]. Note that x < y < z is used instead of x > y > z, refer to (|Cox et al 
19971 Chapter 2, Definition 3 (Lexicographic Order)) for the difference. 



In the case where the polynomials in the input set have many (say more than 6) 
variables, Groebner needs to deal with many intermediate polynomials (see examples 
3, 14, 16, and 17 in Table 2) and CharSet and CharSetw need to deal with polynomials of 
complex initials and high degrees, while NewCharSetw can partially avoid such problems. 
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Finally, we point out that admissible reductions may also be incorporated into other 
algorithms of triangular decomposition, which would create plenty of room for improve- 
ments on all such algorithms. 
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Appendix 

In this appendix we provide more details about the outputs of the five algorithms 
in comparison, including degree tuples, numbers of terms, head monomials, and maxi- 
mal digits of coefficients of the output polynomials, for some (but not all, due to space 
limitation) of the examples. 

Table 3 collects the degree tuples of the polynomials in the output of each algorithm. 
Let us explain the meanings of the entries in this table using an example: the entry in 
the second row and the third column containing three tuples implies that the output of 
CharSet for example 2 consists of three polynomials, say F,G,H with F < G < H, the 
degrees of F in w, x, y, z (w < x < y < z) are 520, 42, 0, respectively. One can see that 
for most examples, the output polynomials of NewCharSetw and NewCharSet are of lower 
degrees than those of CharSetw and CharSet. The entry in the second row and the last 
column contains . . ., which means that there are more elements. 

Statistical data about the numbers of terms of expanded polynomials in the output 
of each algorithm are given in Table 4. The entry 6,6,6 in the second column and the 
second row indicates that for example 1, CharSetw outputs a set of three polynomials, 
which consists of 6, 6, 6 terms respectively. For example 3 there are four polynomials in 
7 variables. The output of any of the first four algorithms contains four polynomials, 
and for NewCharSet some of the polynomials are of degree > 100. This may explain why 
NewCharSet is slower than CharSet. Similar situations occur for examples 8 and 16. For 
examples 2, 6, and 12, the polynomials in the output of NewCharSet have fewer terms 
than those in the output of CharSet, so NewCharSet is faster. For most examples, the 
polynomials contained in the output of NewCharSetw have fewer terms, so NewCharSetw 
is faster than CharSetw, CharSet, and NewCharSet. For examples 3, 14, and 15, Groebner 
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Table 3. Degree tuples of polynomials in the outputs of five algorithms 



No 


CharSetw 


CharSet 


NewCharSetw 


NewCharSet 


Groebner 


Z 




[OZU, 4Z, u, UJ 


To in n nl 
[O, 1Z, U, UJ 


ro 1 9 n nl 

[O, 1Z, U, UJ 


To in n nl 
[O, 1Z, U, UJ 






\a on 1 nl 
[D, OU, 1, UJ 


\a f\ i nl 

[4, O, 1, UJ 


\A a i nl 
[4, 0, 1, UJ 


Ti a 1 1 i nl 

[14, 11, 1, UJ 






[520, 40, 0, 1] 


[2,1,1,1] 


[4, 6, 0, 1] 


[10, 10, 1,0],... 





\ac\ n n nl 

[4U 3 U, U, UJ 


T7Q n n nl 

1 1 O, Uj Uj UJ 


fQfi n n nl 

[OO, u, u, UJ 


[on n n nl 

[OU, u, u, UJ 


\oa n n nl 

[Z4, U, U, UJ 




[oy, i, u, uj 


rero 1 n nl 

[OO, 1, u, UJ 


roc; i n nl 

[ZO, 1, u, UJ 


roe: -i n nl 

[ZO, 1, u, UJ 


roQ i n nl 

[ZO, 1, u, UJ 




[oy, u, ±, uj 


\7f\ n 1 nl 

[ ( U , U , 1 , UJ 


\q 9 i nl 

lO, Z, 1, UJ 


[90 n i nl 
[zy, u, i, uj 


roQ n i nl 

[ZO, U, 1, UJ 




[1,1,1,1] 


[71, 0, 0, 1] 


[1,1,1,1] 


[29, 0, 0, 1] 


[23, 0, 0, 1] 




ri a « n n nl 

[140, U, U, UJ 


ri fifi n n nl 
[IDo, U, U, UJ 


T9a n n nl 
[Zo, U, U, UJ 


\o7 n n nl 
[z < , U, U, UJ 


T9/i n n nl 
[Z4, u, U, UJ 




ri a a i n nl 

[1^, J. j U, UJ 


M A(\ i n nl 

[l^iU, 1 , U, UJ 


[9« 1 n nl 

[ZU, 1, U, UJ 


rofi i n nl 

[ZU, 1 , U, UJ 


roQ i n nl 

[ZO, 1 , U, UJ 




[145, 0, 1, 0] 


[167, 0, 1, 0] 


[1, 3, 1, 0] 


[26, 0, 1, 0] 


[23, 0, 1, 0] 




Iia^ n n 11 

IX^Uj U, U, -LJ 


M fifi n n 11 

[1UU, U, U, 1J 


[0, 2, 0, 1] 


[9fi n n 11 

[ZU, U, U, 1J 


[9^ n n 11 

[ZO, U, U, 1J 






[five; n nl 

[U I o. u, UJ 


ri 1 9 n nl 

[1 1Z, U, UJ 


[Qi n nl 

[Ol , U, UJ 


\ka n nl 

[Ol, U, UJ 






[608, 1, 0] 


[102, 1, 0] 


[80, 1, 0] 


[53, 1, 0] 






\(\7A n il 

1 U I U, 1J 


\A 7 11 


[80, 0, 1] 


r^ n il 

[OO, U, 1J 


10 


ri n nl 

[100, u, UJ 


n n nl 

[lOo, U, UJ 


rifii n nl 

[1U1, u, UJ 


Mfii n nl 

I1D1, u, UJ 


ri a~\ n nl 

[141, U, UJ 




ri 9n i nl 

1 1ZU, 1 , UJ 


ri 9n i nl 

[1ZU, 1 , UJ 


ri9s i nl 

[1ZO, 1 , UJ 


[ioq i nl 

[1ZO, 1, UJ 


ri An 1 nl 

[11U; 1, UJ 




fl 5 11 


[125, 0, 1] 


[5 1 11 

[U, -L, Xj 


[133, 0, 1] 


[140, 0, 1] 


16 


[16,0,0,0] 


[14,0,0,0] 


[16,0,0,0] 


[14,0,0,0] 






[8,1,0,0] 


[8,1,0,0] 


[8,1,0,0] 


[8,1,0,0] 






[2,1,1,0] 


[9,0,1,0] 


[2,1,1,0] 


[9,0,1,0] 






[9,0,0, 1] 


[9, 0, 0, 1] 


[1,1,1,1] 


[9,0,0, 1] 





outputs more polynomials of more terms than NewCharSet. It is therefore slower than 
NewCharSet. 

Table 5 displays the head terms of the polynomials in the output of each algorithm. 
Because most of the initials are quite complicated, we have reproduced only the head 
terms instead of the initials. Table 6 shows the maximal digits of coefficients of the 
polynomials in the output of each algorithm. One sees clearly that the polynomials in 
the outputs of NewCharSet and NewCharSetw have smaller coefficients. 
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Table 4. Numbers of terms of polynomials in the outputs of five algorithms 



No CharSetw CharSet NewCharSetw NewCharSet Groebner 



1 


6,6,6 


6,6,8 


6,6,6 


6,6,8 




2 




5450,93, 10390 


23,12,5 


23, 12, 17 


23,45,35,39,5, 
30,31,41,3 


3 


18,26,38,2 


18, 26, 38, 33 


18,26,5,2 


18,26,114, 108 


10,8,5,17,25,9,... 


1 






3,4,4,6,6,4 


3,4,4,4,4,4 


3,3,3,3,3,3 


5 


41,80,80,4 


74, 107, 142, 143 


37,51,12,4 


37,51,60,60 


25,25,25,25 


6 


143, 284, 284, 286 


165, 288,327,328 


25,50,8,5 


25,50,48,50 


25,25,25,25 


7 


292,522 


292, 522 


88, 158 


88,158 


88, 393, 498 


8 


128,21,35,4 


89,21,35,35 


253,116,7,4 


3217, 116, 152, 15087 


89,220, 119, 150, . . . 


9 






28,55,8, 10,5 




4,8,10,5, 12,13, 
5,12,15, 13,5 


10 


57, 112, 112, 112 


103, 153, 180, 204 


38,137, 45,6 


60, 118, 118, 118 


25,26,27,25,27 


11 




445,792,887 


70, 140,8 


54, 140, 41 


54, 55, 55 


12 




663,1188, 1324 


108, 198, 13 


81, 162, 162 


53,55,55 


13 






3752,29,29,12, 
8,6,8,6,7,6 




172,11,10,11,11, 
10, 11, 11, 11, 10 


11 






87, 199, 16 


76, 152, 150 


45,48,48,48, . . . 


15 


139,153,3 


139, 153, 152 


90, 88, 3 


90,88,88 


90, 142, 142 


16 


3123, 210,8,547 


1984, 210, 378, 547 


20933,210,8,6 


8169,210,378, 547 




17 






49226,305,8,6 
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Table 5. Head terms of polynomials in the outputs of five algorithms 



No 


CharSetw 


CharSet 


NewCharSetw 


NewCharSet 


Groebner 


1 


/, y 2 i,y 2 j 


/, y 2 i,y 3 j 


/, y 2 i, y 2 j 


f,y 2 i, y 3 j 




2 




tu 508 a .42 )a .29 2/j 
w 510 x 40 z 


x 12 ,w 2 x 3 y, 
xyz 


x 12 ,w 2 x 3 y, 
x 6 z 


19 14 
x , w y, 

w 10 xy, . . . 


5 


X 40 , X 39 X2, 
39 

xf*X3, X4 


/v>73 ™,52~. 
1 '1 ^ ' 

70 70 


™36 ™24, r> 
1 '1 -si 

X?,X3, X4 


X 30 ,X 24 X2, 
X 29 X3,X 29 X4, 


X 24 , X2, X3,X4 


6 


d 146 ,d 144 p, 
d 145 c,d 145 q 


d 166 cd 166 (J 


d 2S ,d 26 p, 
p 3 c, pg 


d 27 ,d 2e p, 

d 2(i c, d 26 q 


d 24 ,p,c,q 


8 


t 9 d 2 ,t 5 z, 
t 6 y,tx 


i 9 d, i 5 2, 
t 6 y,t 6 x 


t 12 da 2 ,t 6 daz, 
zy,yx 


t 10 a 3 b 2 c 6 d 5 ,t (i adz, 
t 7 ady,t 8 a 3 b 2 c 6 d 11 x 


t s d,c 2 d 4 z, 
tc 2 d 2 z,t 2 bc 3 z, . . . 


10 


q 56 ,q 55 c, 
q 55 p, q 55 d 


<7 104 ,9 78 c, 

g 91 P ,g 103 rf 


q 73 ,q 62 c, 
q 4 c 2 p, qpd 


q 60 ,q 59 c, 
q 59 p,q 59 d 


24 

c 3 ,p,d 


12 




x 675 ,x 608 t/, 
a; 674 2 


x 112 ,x 102 y, 
x 3 y 5 z 


x sl ,x so y, 
x 80 z 


54 

x ,y, z 


15 


x 153 ,x 120 y, z 


x 15s ,x 120 y, 
x 125 z 


x 161 ,x 12s y, 
x 5 z 


x 161 ,x 12s y, 
x 133 z 


141 


16 


^71 1 ^71 *^72 ; 


14 7 
"^71 ' ■^71*^72 ; 


■^71) -^7x ^72 ? 

3^2:73, a^io^aigo 


14 7 

X 7 J,X7 1 X72, 

^71^73, a; 7 i^80 




17 






xj 6 , xjx2, 
X 2 X3, X1X3X4 
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Table 6. Maximal numbers of digits of coefficients of polynomials in the outputs of five algo- 
rithms 



No 


CharSetw 


CharSet 


NewCharSetw 


NewCharSet 


Groebner 


1 


1,1,1 


1, 1, 1 


1,1,1 


1,1,1 




2 




212,3,212 


2,1,1 


2,1,1 


2,3,3, .. . 


3 


2,2,2,1 


2,2,2,2 


2,2,1,1 


2,2,4,4 


1,1,1,... 


1 






5,15,5,7,7,3 


5, 15, 14, 24, 23, 13 


5,5,4,5,4,3 


5 


40,110, 146,1 


71,55,70,70 


20,27,4,1 


20,28,215,216 


20, 257, 258, 245 


6 


158, 164,776,468 


184, 164, 188, 189 


39,36,2,1 


18, 22, 614, 298 


18,266,271,266 


7 


11, 11 


11, 11 


6,7 


6,6 


6, 15, 17 


8 


2,2,2,1 


2,2,2,2 


5,4,1,1 


6,2,2,9 


2,3,3, .. . 


9 






14,21,1,1,1 




3,3,6, .. . 


10 


36, 123, 330, 383 


67, 50, 60, 288 


20, 160, 9, 1 


72, 264, 603, 904 


15, 163, 186, 143, 182 


11 




94,85,97 


18, 17, 1 


8,15,298 


8, 202, 202 


12 




89,78,99 


16, 14, 1 


12, 22, 364 


9, 193, 194 


13 






3,1,1,1,1,1,... 




2,1,1,1,1,1,... 


14 






17,19,1 


16, 64, 686 


10, 174, 170, . . . 


15 


13,11,1 


13, 11, 11 


9,7, 1 


9,7,7 


9, 625, 626 


16 


6,3,1,3 


5,3,3,3 


8,3,1,1 


6,3,3,3 




17 






8,3,1,1 






18 










6,17, 19,... 
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